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' Gapless superconductivity can arise when pairing occurs between fermion 

(-H , species with different Fermi surface sizes, provided there is a sufficiently large 

^ ■ mismatch between Fermi surfaces and/or at sufficiently large coupling con- 

q ' stant. In gapless states, secondary Fermi surfaces appear where quasiparticle 

excitation energy vanishes. This work focuses on homogeneous and isotropic 
superfiuids in the s-wave channel, with either zero (conventional superconduc- 
tor), one, or two spherical Fermi surfaces. The stability conditions for these 
candidate phases are analyzed. It is found that gapless states with one Fermi 
surface are stable in the BEC region, while gapless states with two Fermi sur- 
faces are unstable in all parameter space. The results can be applied to ultracold 
C ~_ 3 ■ fermionic atom systems. 
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1. Introduction 



This work focuses on stability conditions and their possible relations to 
different Fermi surface topologies in a superconductor with unequal num- 
ber densities of fermions (or unequal chemical potentials) 1 . The physics of 
q | paired fermion systems with unequal densities of the two fermion species 

is of interest in the study of two physical systems: (1) Ultracold atomic 
fermionic gases, where one can freely choose populations in two hyperfine 
states of the fermionic atom. Experimental work is currently being con- 
■ ducted with these systems. (2) Quark matter in the interior of neutron 

stars, which is believed to be a color superconductor. There, the mismatch 
in quark Fermi surfaces is driven by differences in quark masses and electric 
and color neutrality conditions. In both cases, the ground state of such a 
system is the subject of current debate. The study of fermionic superfiu- 
ids with imbalanced populations is a novel subject, because until recently, 
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experimental evidence for this situation was elusive. In a superconductor, 
an imbalance in population between spin-up and spin-down electrons can 
be created by a magnetic field. However, due to the Meissner effect, the 
magnetic field is either completely shielded from the superconductor bulk 
or enters in the form of quantized flux lines or vortices. 

Cold fcrmionic atomic gases provide new possibilities in the experimen- 
tal exploration of fermionic superfluids with unequal mixtures of fermions. 
In cold atoms, Feshbach resonance permits tuning of the interaction from 
(I) attraction and a resulting superfluid of loosely bound pairs (the BCS 
regime at £ > n^ 1 / 3 ); to (2) repulsion and a resulting Bose-Einstein con- 
densate of tightly bound molecules (the BEC regime at C < n -1 / 3 ), where £ 
is the size of a pair and n~ 1//3 is the interparticle separation (or, more pre- 
cisely, the mean free path). With equal mixtures, the BEC-BCS crossover is 
smooth, with no phase transition. With asymmetric densities, one or more 
phase transitions are expected, and a more complex phase diagram may 
result. 

In the weak coupling regime, a BCS superfluid remains stable as long as 
the difference in chemical potential is small compared to the pairing gap, 
Sfi < A; the gap prevents the excess unpaired atoms from entering the 
superfluid state. By either increasing the mismatch or reducing the binding 
energy (and hence decreasing the gap), a quantum phase transition from the 
superfluid to normal state takes place, and superflidity ceases. The point at 
which the phase transition occurs is known as the Clogston limit, which can 
be estimated as Sfi ~ A ~ ^exp(— 1/g) <C fi, where (i is the Fermi energy. 
Thus only an exponentially small population imbalance is allowed in weak 
coupling. In strong coupling, as one approaches Feshbach resonance, the 
situation is quite different. On the BEC side, superfluidity with imbalanced 
population is robust over a wide range of parameter space. Surprisingly, on 
the BCS side, macroscopic imbalance is also possible in a superfluid state, 
possibly due to the formation of a gapless superfluid at strong coupling. 
This gapless superfluid incorporates large numbers of unpaired fermions 
which reside at the secondary Fermi surfaces. Here, we consider stability 
conditions for gapless states with different Fermi surface topologies. 

2. Definitions for the screening masses and susceptibility. 

We consider two species of nonrelativistic fermions, V — C^iiV^) with 
the same mass but with different chemical potentials, fi — diag(^i, /i 2 )- 
There is an attractive interaction only between different species of fermions, 
g (^(72^*) i^ T ^2' l f)- The order parameter is & = Aa 2 where A = 2g < 
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i\) T (ji^ >, it defines pairing in the singlet channel, spin 0, [2] x [2] = [l] + [3]. 
The order parameter breaks the original group associated with conservation 
of particle numbers of each species down to the diagonal subgroup associ- 
ated with conservation of difference in particle numbers, U(l) ai x U(l) a2 — > 
U(l) ai - a2 which is invariant under simultaneous rotation a\ = —0:2- To- 
tal number of particles is not conserved. Formally, one can express this 
pattern of symmetry breaking by gauging the theory. We introduce two ex- 
ternal gauge fields, and couple each species of fermions to its gauge field, ipi 
couples to A\ with g\ and -02 couples to A 2 with g 2 , that is reflected by the 
generators of the gauge group, T\ = diag(l, 0) and T 2 — diag(0, 1). Accord- 
ing to the symmetry breaking pattern, gauge fields mix in a superconductor, 
rotated sum of fields A\ = A1+A2 is screened with non-zero Meissner mass 
via Andersen-Higgs mechanism, m 2 M 7^ corresponding to U(l) ai + a3 , and 
rotated difference of fields A 2 = A\ — A 2 propagates in a superconductor 
unscreened, Meissner mass is zero m 2 M = corresponding to U(l) ai - a2 . 
There is a different mixing of gauge fields in case of Debye masses. In QCD, 
with color and electromagnetic gauge groups, SU(3) C x U(1)em, mixing of 
gauge fields is more complicated. 

We treat four-fermion interaction on a mean field level. We integrate 
over the fermions, obtain fermion determinant, 



VA exp 



IAI 2 1 
S A + ^---Trln(S- 1 +A) 



(1) 



We use the Nambu-Gorkov formalism with particle-hole basis, * = {ip, ip*). 
Inverse fermion propagator and gauge field are 4x4 matrices in this basis, 
including the Nambu-Gorkov and two fermion species indices, 



[GoT 1 
* + [Go] 



(2) 



where the inverse free fermion propagators are [Gq] 1 = idt ± ^ ± /U, 
and we have abbreviated for the gauge fields A = diag(A + ,A~) with 
A ± = ±T a A° a T £ - ^(V • A a + A a • V) and T a = g a T a . Perform- 
ing derivative expansion and collecting terms quadratic in the gauge field, 
we produce fermion loops, II™, n^, 11^, W^ b . Debye mass in one- loop 
is defined by the temporal component of polarization operator, Meissner 
mass is given by the spatial component, m 2 D ab = — lim p ^o n™(0, p) and 
m M,ab = \ nm P ^o(% - PiPj)n^ 6 (0, p), where p t = pi/p. Screening masses 
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are given 1 



™l,ab =-|m~X; Tr[S(K)T-S(K P)T^\ , 

~^ K 

+ £ [1 ~ {p ' k)2] T ^ s ( K ) T aS(K P)r+]) , 

where S 1 is the fermion propagator and we have introduced the follow- 
ing matrices in Nambu-Gorkov space, = diag(T a ,±r a ) and t 2 a = 
diag(r^, — r„). Physically, Debye mass to all loops, including ladder dia- 
grams, is equivalent to compressibility of a system, and Meissner mass can 
be associated with the density of superconducting fermions. These 2x2 
matrices in two-fermion space shall be evaluated in the following section in 
order to obtain stability conditions for gapless superconductors. 

One may derive the pressure from the partition function in (1) using 
the Cornwall- Jackiw-Tomboulis formalism. The pressure is the negative 
of the effective potential at its stationary point (i.e., with the propaga- 
tors determined to extremize the effective potential). The fcrmionic part 
of the pressure is p = i^Trln^" 1 + i^Tr^S - 1] + T 2 [S], where 
So = diag(Gp ,Gq ) is the tree- level fermion propagator in Nambu-Gorkov 
space and ^[5] is the sum of all two-particle irreducible diagrams. The 
number densities is defined n a = dp/dfx a . The number susceptibility \ is 
defined as the derivative of the number density with respect to the chemical 
potential (at constant volume and temperature). Using the expression for 
the pressure, we obtain 1 



1 T 



(4) 



where £ is the fermion self-energy, S^ 1 = Sq 1 + X. The first term on 
the right-hand side of this equation is given by the one-loop result for the 
electric screening mass, cf. Eq. (3). For the second term, we assume that 
the self-energy £ depends on fj, only through the gap, then we obtain 1 

m2 p,ab . dn a dA 

Xab = h ~jrrj; — ■ (5) 

g a g b dA d^ b 
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In general, the self-energy X contains terms of any number of fcrmion loops. 
Consequently, the number susceptibility contains terms of arbitrary many 
fermion loops too, corresponding to the exact Debye mass including all 
possible perturbative insertions. Remarkably, the free fermion result for \, 
i.e. £ = 0, gives the one- loop result for m 2 D . Susceptibility is equivalent to 
compressibility of a system. We shall use Eq. (5) in the following section to 
compute the number susceptibility. As this equation shows, it goes beyond 
the one-loop result for the electric screening mass. 

3. Results 

We consider three cases distinguished by how many zeros quasiparticle dis- 
persion has, = \J (k 2 /2m — p) 2 + A 2 — Sfj,, where /x is the average chem- 
ical potential and 6/j, is the difference between potentials, Fig. 1. In this 
case, number of zeros is equivalent to number of the Fermi surfaces. No 
zeros corresponds to the fully gapped state. In case of one Fermi surface, 
momenta outside the Fermi surface contribute to the pairing, while the ex- 
cess of fcrmions resides inside the Fermi ball. In case of two Fermi surfaces, 
the excess of fcrmions resides between the two Fermi surfaces in momentum 
space. As was noticed by Son, apart from number of zeros, dispersion can 
have two different characteristic behaviors, distinguished by the position 
of the minimum. Minimum is located at nonzero momentum for positive 
p, and corresponds to BCS, and minimum shifts to p = for negative fi, 
and corresponds to BEC. This behavior manifests itself in stability condi- 
tions. We depict different topologies which are distinguished by the number 
of Fermi surfaces on the phase diagram in dimcnsionlcss average chemical 
potential p/A and difference in chemical potentials 5fi/A, the gap is the 
energy scale, in Fig. 2 regions between the solid lines. We have fully gapped 
F 0} and gapless states with one F\ and two F 2 Fermi surfaces. At small mis- 
match, there is a fully gapped state F , which is at positive p and Sfi < A, 
BCS, and at small and negative p, p < — y/S/2 2 — A 2 , BEC. Increasing 
mismatch, when mismatch is at least larger than the gap, there is a gapless 
state with one Fermi surface F\ when fj, is around the Feshbach resonance, 
— \J 8[i 2 — A 2 < p < \J S[i 2 — A 2 , and gapless state with two Fermi surfaces 
F2 for positive fj, restricted from below, y/ S/i 2 — A 2 < p. We therefore ex- 
pect that Fi exists at strong coupling, while F 2 probably exists only at 
weak coupling. 

The question we are solving here is, "What is the ground state in a 
degenerate Fermi system with asymmetric number densities of fermions?" . 
We avoid solving for the ground state explicitly. Instead, we check stabil- 
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Fig. 1. Schematic plot of possible quasiparticle occupation numbers (in arbitrary units). 



ity criteria, positive definite eigenvalues of screening masses and number 
susceptibility. If stability conditions are satisfied, homogeneous supercon- 
ducting gapless state is indeed the ground state in this parameter range. If 
stability conditions are not satisfied, the alternative state may be realized 
as the ground state. These include LOFF, spatially separated mixtures, 
normal (not superconducting) states. 

We analyze stability conditions, Eq. (3) and Eq. (5), in all parameter 
space (/2,5/i), and depict stable/unstable regions together with topology 
regions, Fig. 2 left panel. Debye mass e.v. are positive in all parameter 
space, hence Debye mass does not impose a constraint. All entries for the 
Meissner mass matrix are the same, it is trivial to diagonalize. We obtain 1 
m M = corresponding to the unbroken sector U(l) ai — ai , and m 2 M = 2L 

3/2 3/2 

corresponding to the broken group U(l) ai + a2 , where L = I + ~ , 

2jj- v /»? 2 -1 

which defines the dashed-dotted (blue online) curve in Fig. 2 left panel. It 
renders all states between the dashed-dotted curve and the solid vertical line 
unstable. There is a strip left in gapless superconductor state with two Fermi 
surfaces F2 which is stable, gray area. Fq and Fx are stable with respect to 
m 2 M everywhere. In magnetic sector, mixing does not depend on chemical 
potentials and it is defined by the pattern of symmetry breaking. In electric 
sector, mixing depends on chemical potentials, e.g. in QCD mixing depends 
on Sp. Mixing in electric and magnetic sectors is the same only for the fully 
gapped case. 

Analyzing the number susceptibility matrix, Eq. (5), we obtain 1 the 
expression defining the sign of e.v., which is very similar in structure to 

I/ 2 , 1/2 

that for the Meissner mass, see expression for L, R = I - p+ , p - . Here, 

^vV- 1 

we defined I = I p (0, 00) — I p (^/pZ, y/p+) and I = I P (0, 00) — I p (y/pZ, yfp+) 
with elliptic integrals I p (a,b) = J dx x 2 /[{x 2 — p) 2 + l] 3 / 2 and I p (a.b) = 
J^dxx 4 " /[(x 2 ~ p) 2 + I] 3 / 2 , and p± = p ± \] r\ 2 - 1 are zeros of = 
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with p = /i/A, rj = 5/i/A. In expression for R one should put p_ = 
when applied to F\ region. It defines the dashed (red online) line, and it 
renders all states between dashed and solid vertical lines unstable. Thus all 
Fi states are unstable, and there is a strip in F\ which is stable, dark gray 
area. F is stable everywhere. Stable states with respect to \ correspond to 
local maximum of pressure. 

We obtained stable regions which are local maxima of pressure, Fig. 2 
left panel. Now we consider global maxima, Fig. 2 right panel. For this we 
compare pressure of the superconducting and normal states. Superconduct- 
ing states which pressure is higher than that of the normal state are stable, 
Ap = p s — p n > 0. The dashed-dotted (blue on line) line is Ap = 0, it 
renders all states above and to the right of it unstable. All unstable regions 
with respect to the screening masses and number susceptibility are subset 
of unstable region with respect to the pressure. In a weak coupling, the ver- 
tical dashed-dotted line reproduces the known Clogston limit, S/u, = 
above which BCS is unstable. The global stability line cuts through the 
stable strip of F\ state, below is a stable superconducting state, above is a 
metastable state. Both lines coincide at large mismatches. 

Currently, experiments are being performed with unequal mixtures of 
fcrmions to map superfluid regions as a function of population imbalance, 
interaction strength and temperature. The experimental signature of su- 
perfluidity is the existence of vortices, which prove phase coherence in a 
sample. 
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